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ABSTRACT 

A  series  of  monthly  indices  of  intensity  of  large-scale,  wind-induced  coastal  upwelling  at 
selected  locations  along  the  west  coast  of  North  America  is  presented  for  the  period  1946 
through  1971.  The  indices  are  based  on  calculations  of  offshore  Ekman  surface  wind  transport 
from  monthly  mean  surface  atmospheric  pressure  data.  Summaries  by  quarter  and  by  year  are 
included. 

The  effect  of  using  monthly  mean  pressure  data  is  evaluated  by  comparison  to  a  similar 
series  of  monthly  means  of  transports  computed  6-hourly  over  a  43^-yr  period.  The  correla- 
tion between  the  two  series  at  an\  particular  location  was  found  to  be  high.  However,  the  slope 
of  the  regression  line  varies  at  tlifferent  locations.  Also  values  off  southern  California  during 
summer  may  be  amplified  relati\  e  to  other  locations  as  a  result  of  the  discontinuity  in  the  atmos- 
pheric pressure  gradient  caused  by  the  coastal  mountain  range  between  the  thermal  low  in 
the  interior  of  southern  California  and  the  higher  pressure  offshore.  The  conclusion  is  that  these 
series  may  be  satisfactory  indicators  of  temporal  variations  of  upwelling  at  each  location,  but 
less  satisfactory  indicators  of  spatial  distributions. 


INTRODUCTION 

Coastal  upwelling  can  profoundly  affect  the 
physical  environment  of  marine  organisms.  An 
extreme  example  is  the  "El  Niiio"  of  the  Peru 
Current  where  high  surface  temperatures  re- 
lated to  cessation  of  upwelling  have  apparently 
resulted  in  such  dislocation  of  marine  fauna  as 
to  have  caused  starvation  of  great  numbers  of 
guano-producing  seabirds  (Wooster,  1960).  Per- 
haps even  more  important  from  a  fishery  stand- 
point is  the  role  of  upwelling  in  the  cycle  of  organ- 
ic production  in  the  sea,  both  in  transporting 
nutrient-rich  deeper  waters  into  the  illuminated 
surface  layers  where  they  are  available  for  pro- 
duction, and  in  providing  virgin  water  which  is 
sufficiently  free  from  predators  to  allow  accumu- 
lation of  large  phytoplankton  blooms  (Gushing, 
1969).  Estimates  have  run  as  high  as  Ryther's 
(1969)   suggestion   that   upwelling   areas   com- 


'    Pacific   Environmental  Group.    National   Marine   Fish- 
eries Service.  NOAA.  Monterey.  CA  93940. 


prising  about  one-tenth  of  1%  of  the  ocean  surface 
may  produce  one-half  of  the  world's  harvestable 
fish  supply.  Evidently,  some  indication  of  fluctua- 
tions in  the  upwelling  regime  is  extremely  impor- 
tant to  understanding  fluctuations  in  marine  pop- 
ulations. This  report  is  an  attempt  to  provide 
indices  of  intensity  of  coastal  upwelling  off  west- 
ern North  America  on  time  and  space  scales  use- 
ful to  fishery  research  in  general,  and  to  the 
MARMAP  (Marine  Resources  Monitoring,  As- 
sessment and  Prediction)  program  of  the  Na- 
tional Marine  Fisheries  Service  in  particular. 

METHOD 

It  is  generally  recognized  that  coastal  upwell- 
ing in  broad,  diffuse  eastern  boundary  currents 
such  as  the  California  Current  and  the  eastern 
portion  of  the  Gulf  of  Alaska  Gyre  is  largely  due 
to  replacement  from  below  of  surface  water 
transported  offshore  by  the  stress  of  the  wind  on 
the  sea  surface  (Sverdrup  et  al.,  1942,  p.  501). 
Other  possible  mechanisms  producing  coastal  up- 


welling  include  upwarping  of  density  surfaces 
within  a  geostrophic  current  (Stommel  and 
Wooster,  1965)  and  cross-isobar  onshore  flow  at 
the  bottom  due  to  bottom  friction  (Hsueh  and 
O'Brien,  1971).  However,  these  are  more  likely 
to  be  important  in  narrower,  more  intense  flows 
than  normally  occur  in  the  area  covered  by  this 
report. 

Our  present  understanding  of  surface  wind 
transport  is  based  on  Ekman's  (1905)  theory. 
Under  Ekman's  assumptions  of  steady  state 
motion,  uniform  wind,  and  infinite  homogeneous 
ocean,  the  mass  transport  per  unit  width  of  ocean 
surface  is  directed  90  degrees  to  the  right  (in  the 
Northern  Hemisphere)  of  the  direction  toward 
which  the  wind  is  blowing  and  is  related  to  the 
magnitude  of  the  wind  stress  by  the  expression 


M 


/ 


(1) 


where  M  is  the  mass  transport  resulting  from  a 
wind  stress,  r,  and  /  is  the  Coriolis  parameter. 
This  mass  transport  has  come  to  be  called  the 
Ekman  transport.  The  layer  in  which  appreciable 
transport  occurs  is  often  referred  to  as  the  Ekmun 
layer  and  extends  from  the  surface  to  depths  not 
exceeding  50  to  100  m.  The  bottom  of  the  Ekman 
layer  is  sometimes  identified  with  the  bottom  of 
the  homogeneous  wind-mixed  zone. 

Smith  (1967)  has  shown  that  Yoshida's  (1955) 
expression  for  the  offshore  transport  in  the  early 
stage  of  coastal  upwelling  reduces  to  the  Ekman 
transport  expression  if  the  stress  is  assumed 
constant.  This  leads  to  the  conclusion  that  the 
Ekman  theory  gives  a  valid  description  of  wind- 
driven  offshore  flow  in  the  early  nonequilibrium 
phase  of  upwelling  as  well  as  in  the  later  steady- 
state  phase. 

The  approach  taken  in  generating  the  indices 
presented  in  this  report  has  been  to  estimate  the 
monthly  mean  wind  stress  on  the  sea  surface  at 
points  near  the  coast,  from  this  to  compute  the 
Ekman  transport,  and  finally  to  resolve  the  com- 
ponent of  Ekman  transport  perpendicular  to  the 
coast.  The  magnitude  of  the  offshore  component 
is  considered  an  indication  of  the  amount  of  water 
upwelled  through  the  bottom  of  the  Ekman  layer 
to  replace  that  driven  offshore  (Fig.  1).  Negative 
values  indicate  onshore  transport  or  convergence 
at  the  coast.  Since  accumulation  of  surface  wa- 
ters tends  to  cause  downward  displacement  of 
the  density  structure  in  the  coastal  area,  this  sit- 


uation is  sometimes  referred  to  as  downirelling. 

The  basic  input  data  is  the  wind  field  over  the 
ocean.  However,  the  distribution  of  sea-surface 
wind  observations  in  the  near  coastal  regions  of 
the  northeastern  Pacific  is  uneven,  both  spatially 
and  temporally.  The  number  available  for  a  given 
area  during  a  given  month  is  often  insufficient  to 
form  a  good  estimate  of  the  monthly  mean  stress 
on  the  sea  surface.  In  order  to  construct  a  consis- 
tent series,  use  is  made  of  the  relationship  in  mid- 
latitude  regions  of  wind  to  atmospheric  pressure. 
Incorporating  atmospheric  pressure  data  in- 
creases the  coverage  in  data-sparse  areas  and 
allows  an  understanding  of  the  behavior  of  large- 
scale  pressure  systems  to  aid  in  providing  contin- 
uity to  the  analysis  of  scattered  observations. 
Therefore,  winds  derived  from  analyzed  atmo- 
spheric pressure  fields  are  used  in  the  produc- 
tion of  these  indices. 

Calculations 

The  computed  values  are  based  on  monthly 
mean  pressure  fields  prepared  by  Fleet  Numeri- 
cal Weather  Central  (FNWC).  These  data  are 
available  on  a  63  by  63  point  square  grid  which  is 
superimposed  on  a  polar  stereographic  projection 
of  the  Northern  Hemisphere  (Hughes,  1966).  The 
mesh  length  is  200  nautical  miles  at  lat.  60  °N  and 
decreases  southward  to  about  144  nautical  miles 
at  lat.  20  °N.  The  data  were  transferred  to  a 
3-degree  mesh  length  geographical  (spherical 
coordinates)  grid  (Fig.  2)  using  Bessel's  central 
difference  formula. 

First  derivatives  of  the  surface  pressure  at 
each  grid  point  were  estimated  by  taking  the  dif- 
ference in  pressure  between  the  grid  points  to 
either  side  and  dividing  by  the  6-degree  angular 
mesh  length.  For  example,  the  derivatives  of  the 
pressure  at  point  "0"  in  Figure  3(a)  would  be  ap- 
proximated as 


dP_ 
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2h 
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where  6  and  A  denote  the  northward  and  east- 
ward angular  coordinates,  h  is  the  3-degree  an- 
gular mesh  length  in  radians  and  P,  denotes  the 
pressure  at  point  "1",  etc.  The  geostrophic  wind 
was  computed  according  to 

I        dP  I  dP 
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Figure  1.  — A  conceptual  diagram  of  the  coastal  upwelliriK  process.  The  coast  of  the 
continent  is  represented  in  cutaway  view  with  the  ocean  to  the  left  of  the  figure.  Off- 
shore transport  in  the  surface  Ekman  layer  due  to  stress  of  the  wind  parallel  to  the 
coast  on  the  sea  surface  is  replaced  by  upwelling  from  depth. 


where  Vg  and  Uy  are  the  respective  northward 
and  eastward  components  of  geostrophic  wind 
velocity,  /  is  the  Coriolis  parameter,  p„  is  the 
density  of  air  and  R  is  the  mean  radius  of  the 
earth.  The  density  p„  was  considered  constant  at 
0.00122  g  cm  ' .  An  estimate  of  the  wind  near  the 
sea  surface  was  formed  by  rotating  the  geo- 
strophic  wind  vector  15  degrees  to  the  left  and 
reducing  it  by  30%  as  in  Figure  3(b)  to  approxi- 
mate frictional  effects. 

The  sea-surface  stress  was  then  computed  ac- 
cording to  the  classical  square-law  formula 


(4) 


where  t  is  the  stress  vector,  p,;  is  the  density  of 
air,  Cfi  is  an  empirical  drag  coefficient,  v  is  the 
estimated  wind  vector  near  the  sea  surface  with 
magnitude  \r\.  A  relatively  high  value,  0.0026, 
of  the  drag  coefficient  was  used  to  partially  off- 
set the  effect  of  using  mean  data.  Finally,  the 
Ekman  transport  was  computed  according  to 
Equation  (1).  Except  for  the  estimates  of  the 
pressure  derivatives  which  necessarily  differ 
because  of  the  different  grid  format,  this  calcu- 
lation procedure  is  the  same  as  that  used  by 
Fofonoff  (I960).-  The  offshore  transport  is  de- 
termined by  resolving  the  component  perpendic- 


ular to  a  line  drawn  by  visual  estimation  on  a 
bathymetric  chart  along  the  dominant  trend  of  a 
200-mile  segment  of  coastline  centered  near  the 
grid  point  in  question. 

The  indices  generated  by  this  method  must  be 
considered  as  indicative  of  rather  large-scale 
coastal  upwelling.  The  6-degree  gap  across  which 
the  derivatives  are  measured  as  well  as  the  scale 
of  coastline  resolution  indicates  a  100-  to  200-mile 
scale  as  being  appropriate  for  interpretation  of 
the  indices.  Certainly,  smaller  scale  upwellings 
could  occur  within  a  larger  scale  downwelling 
just  as  short-term  upwellings  might  occur  during 
a  monthly  period  of  average  downwelling. 

Pressure  Data 

The  number  of  synoptic  surface  atmospheric 
pressure  fields  available  at  FNWC  for  each 
month  during  the  period  1946-71  is  variable. 
Because  of  the  nonlinear  dependence  of  stress  on 
wind  velocity  (Equation  4),  use  of  a  variable  sam- 
pling interval  destroys  the  internal  consistency 
of  the  time-series.  For  this  reason  monthly  mean 


-  Fofonoff.  N.  P.  1960.  Transport  computations  for  the 
North  Pacific  Ocean-1958.  Fish  Res.  Board  Can.,  Manuscr. 
Rep.  Ser.  (Oceanogr.  Limnol.)  80,  81  p.  (Processed.) 
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Figure  2.  — Data  grid.  Intersections  at  vvhicli  upwelling  indices  are  computed  are  marked  with  large  dots. 


♦  * 


.0*  ,^ 


(a) 


07- 


GEOSTROPHICt  WIND 


(b) 


Figure  3. —  (a)  Configuration  of  data  points  used  to 
form  pressure  derivatives,  lb)  Diagram  showing 
the  transformation  of  the  geostrophic  wind  vector 
to  form  an  estimate  of  the  wind  near  the  sea  surface. 


pressure  fields  were  used  as  input  data,  thereby 
making  best  use  of  available  data  within  each 
month  while  retaining  a  consistent  sampling  in- 
terval for  computation  of  the  stress.  The  effect  of 
using  monthly  mean  data  on  the  absolute  magni- 
tude of  the  index  at  each  location  is  discussed  in 
the  next  section. 

The  monthly  mean  fields  have  been  assembled 
by  FNWC  from,  as  far  as  was  possible,  6-hourly 
(four  per  day)  synoptic  fields.  Where  this  fre- 
quency was  not  available,  twice  per  day  where 
available  or  once  per  day  fields  were  utilized. 
This  decreased  frequency  occurs  mainly  in  those 
years  previous  to  1963  when  FNWC  began  to 
produce  its  own  analysis.  To  extend  the  series 
back  to  1946,  FNWC  has  collected  synoptic  fields 
from  the  sources  listed  in  Table  1. 

Because  of  the  different  data  sources  and  the 
expected  trend  toward  increased  detail  in  later 
pressure  fields  due  to  an  expanded  system  of 
meteorological  observations,  care  should  be  used 
in  attempting  to  discern  long-term  trends  in 
these  upwelling  indices.  For  example,  it  is  diffi- 
cult to  determine  how  much  of  the  decreased 
amplitude  of  the  yearly  cycle  apparent  at  several 
of  the  locations  during  the  earliest  several  years 
is  real  and  how  much  is  due  to  relative  smoothing 
of  the  data  in  the  earliest  pressure  fields.  The  in- 
tention has  been  to  produce  a  series  which  will 
give  a  useful  indication  of  the  relative  intensitv 


of  upwelling  by  month  within  any  group  of  fairly 
contemporaneous  years. 

Effects  of  Monthly  Mean  Data 
on    Calculated    Results 

Due  to  the  nonlinear  linkage  of  wind  to  stress 
(Equation  4)  the  value  of  the  Ekman  transport 
computed  from  monthly  mean  data  will  be  small- 
Table  1.  — Sources  of  synoptic  surface  atmospheric  pres- 
sure fields  used  by  Fleet  Numerical  Weather  Central  in 
constructing  monthly  mean  fields. 


Time 
period 


Source 
agency 


Remarl^s 


Jan.    iy46-   National  Climatic  Center    Hand  analysis  of  available 
Mar.  1955    Ashville.  N.C.  reports:  once  per  day. 


Apr.  1955-  National  Center  for 
Dec.  1959    Atmospheric  Research 
Boulder,  Colo. 


Jan.    1960-  National  Climatic  Center 
June  1962    Asheville,  N.C. 


National  Meteorological 
Center  analysis  digitized 
at  alternate  gridpoints; 
twice  per  day. 

Once  per  day. 


July    1962-  Fleet  Numerical  Weather    FNWC  objective  machine 
Dec.  1971    Central.  Monterey,  Calif,    analysis.  Certain  gaps 

filled  with  NCC  data; 

four  times  per  day. 


er  than  the  monthly  mean  of  Ekman  transport 
computed  at  synoptic  intervals  by  an  amount 
depending  upon  the  variance  of  the  wind  vector 
within  the  month.  In  order  to  evaluate  this 
effect,  a  54-mo  series  (January  1967  through  June 
1971)  of  monthly  means  of  6-hourly  computations 
of  offshore  component  of  Ekman  transport  was 
generated  for  comparison  with  the  corresponding 
values  computed  from  monthly  mean  data.  The 
results  of  the  comparison  are  summarized  in 
Table  2. 

The  high  values  of  the  correlation  coefficients 
indicate  quite  linear  relationships  between  the 
two  sets  of  series.  Thus  the  distortion  introduced 
into  the  series  at  any  particular  location  due  to 
use  of  monthly  mean  data  is  mainly  one  of  abso- 
lute magnitude.  The  relationship  of  each  monthly 
value  to  the  other  monthly  values  at  the  location 
is  affected  only  slightly.  This  indicates  that  the 
series  computed  from  monthly  mean  data  gives 
an  indication  of  the  monthly  variations  in  inten- 
sity of  upwelling  nearly  equivalent  to  that  given 
by  the  corresponding  series  of  monthly  means  of 
values  computed  each  6  hr. 

However,  Table  2  shows  the  slope  of  the  re- 
gression line  for  the  two  series  to  be  variable 
from  one  location  to  the  other.  This  indicates  a 
change  with  location  in  the  amount  of  variability 
of  the  wind  vector  within  a  month.  Thus  the  use 
of  monthly  mean  data  in  conjunction  with  the 
nonlinear  stress  law  (Equation  4)  gives  an  erro- 
neous impression  of  the  relative  magnitudes  of 
offshore  Ekman  drift  when  different  locations 
are  compared.  For  this  reason  the  indices  tabula- 
ted in  this  report  which  are  designed  to  indicate 
temporal  variations  at  a  particular  location, 
should  be  used  with  caution  as  indicators  of  spa- 
tial   distributions. 

This  distortion  of  the  spatial  distributions  in- 
troduced by  designing  the  computational  method 
to  minimize  distortion  of  the  time-series  at  each 
particular  location  can  be  illustrated  by  exam- 
ining the  long-term  mean  annual  cycle  of  off- 
shore Ekman  transport  computed  in  two  differ- 
ent ways.  Figure  4(a)  is  a  time-series  isogram  of 
offshore  component  of  Ekman  transport  for  a 
long-term  composite  year  computed  from  actual 
wind  observations  taken  from  the  National  Cli- 
matic Center's  file  of  marine  surface  observa- 
tions. The  figure  summarizes  over  75,000  indi- 
vidual wind  observations  taken  by  ships  at  sea 


Table  2.— Results  of  comparison  of  a  54-mo  series  of 
monthly  means  of  offshore  component  of  Ekman  transport 
computed  6-hourly  to  the  corresponding  series  computed 
from  monthly  mean  data  at  the  selected  coastal  grid- 
points. 

Rank  Product-  ^  ^^'^P^^''^  ^ 

Location  correlation    moment  correlation  " 

coefficient  coefficient  ^ 

line 


60°N.149°W 

0.99 

0.98 

0.62 

60°N.146°W 

0.99 

0.98 

0.66 

57°N.137°W 

0.99 

0.98 

0.66 

53°N.134°W 

0.97 

0.94 

0.54 

51°N,131°W 

0.96 

0.94 

0.47 

48°N.125°W 

0.99 

0.97 

0.48 

45°N,125°W 

0.99 

0.97 

0.49 

42°N,125°W 

0.99 

0.96 

0.57 

39°N.125°W 

0.99 

0.96 

0.68 

36°N,122°W 

0.99 

0.98 

0.83 

33°N.119°W 

0.99 

0.99 

0.94 

30°N.1I9°W 

0.99 

0.99 

0.98 

27°N.116°W 

0.95 

0.96 

0.97 

24°N.113°W 

0.96 

0.96 

0.83 

21°N.107°W 

0.99 

0.96 

0.59 

over  the  20-yr  period,  1948  through  1967,  within 
the  1-degree  squares  shown  in  the  accompanying 
coastline  plot.  The  Ekman  transport  was  comput- 
ed for  each  wind  observation  and  the  offshore 
components  were  averaged  by  month  for  each 
coastal  square.  The  drag  coefficient  used  with 
these  actual  observed  winds  was  0.0013.  Bakun 
(1971)  demonstrated  that  the  spatial  distribu- 
tions obtained  in  this  manner  agreed  well  with 
the  mean  annual  cycle  of  sea-surface  temperature 
distributions.  Reference  to  the  slopes  of  the  re- 
gression lines  in  Table  2  indicates  that  since  this 
drag  coefficient  is  half  that  used  with  winds  com- 
puted from  monthly  mean  pressure  fields  to  gen- 
erate the  indices  tabulated  in  this  report,  rough 
numerical  equivalence  between  the  two  trans- 
port sets  should  occur  at  about  lat.  45°N. 

Figure  4(b)  is  a  similar  isogram  constructed 
from  long-term  mean  monthly  values  for  the  same 
period  of  the  indices  derived  from  monthly  mean 
pressure  data.  The  locations  at  which  the  indices 
were  calculated  are  marked  on  the  accompanying 
coastline  plot.  Comparison  with  Figure  4(a)  shows 
an  extreme  southward  shift  and  intensification 
of  the  point  of  maximum  indicated  upwelling. 

Figure  4(c)  displays  the  data  from  Figure  4(b) 
after  having  been  adjusted  to  the  slopes  of  the  re- 
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Figure  4.  — (a)  Long-term  mean  annual  cycle  of  offshore  Ekman  transport  coniiiosilcd  hum  actual 
wind  observations  within  the  1-degree  squares  indicated  on  the  coastline  drawing  to  the  right, 
lb)  Long-term  mean  annual  cycle  of  offshore  Ekman  transport  formed  from  up  welling  indices  com- 
puted froin  monthly  mean  atmospheric  pressure  data  at  points  indicated  by  dots  on  the  coastline 
drawing,  icl  Values  from  Figure  4(al  after  having  been  adjusted  to  the  slopes  of  the  regression 
lines  in  Table  2.  Units  are  cubic  meters  per  second  per  meter  of  coastline. 


gression  lines  by  multiplying  each  monthly  value 
by  0.5  and  by  the  reciprocal  of  the  slope  given  in 
Table  2.  The  multiplication  by  0.5  is  necessary  be- 
cause the  value  of  the  drag  coefficient  used  to 
construct  Figure  4(a)  is  half  that  used  for  Figure 
4(b)  whereas  the  slopes  in  Table  2  were  computed 
for  two  series  with  equal  drag  coefficients.  Multi- 
plication by  the  reciprocal  of  the  slopes  then 
serves  to  adjust  the  data  so  that  the  spatial  varia- 
tion in  the  underestimate  of  the  monthly  mean 
wind  stress  due  to  computation  from  monthly 
mean  atmospheric  pressure  data  as  predicted  in 
Table  2  is  removed. 

The  magnitudes  and  spatial  variations  dis- 
played in  Figure  4(c)  are  more  like  those  shown  in 
Figure  4(a)  than  are  those  in  Figure  4(b).  Some 
smoothing  of  gradients  in  Figure  4(b)  and  (c)  re- 
sults from  having  data  points  at  3-degree  intervals 
with  each  point  incorporating  data  across  a  6- 
degree  area.  This  causes  some  loss  in  detail  and 
displacement  of  contours  relative  to  Figure  4(a| 
where  data  points  are  at  1-degree  intervals  and 
incorporate  data  restricted  to  a  1-degree  square 
area.  There  remains  however  a  definite  shift  of 
the  maximum  toward  southern  California  in  Fig- 
ure 4(c)  relative  to  Figure  4(a).  A  possible  cause 
for  this  shift  in  maximum  upwelling  is  described 
in  the  following  section. 


An  Effect  of  Coastal  Topography 

During  the  summer  an  intense  thermal  low  de- 
velops in  the  interior  of  southern  California.  Due 
to  the  3-degree  mesh  length  of  the  computation 

grid,  the  meridional  component  of  the  geostrophic 
wind  at  a  coastal  gridpoint  is  computed  (Equation 
2)  from  a  continuous  constant  pressure  gradient 
between  an  offshore  gridpoint  and  one  on  the  con- 
tinent (dashed  line  in  Fig.  5).  However,  the  south- 
ern California  coastal  mountain  range  causes  a 
discontinuity  in  this  pressure  gradient  such  that 
the  gradient  actually  in  equilibrium  with  the  geo- 
strophic wind  (solid  line  in  Fig.  5)  may  be  less 
than  that  used  in  the  computations.  This  leads  to 
an  overestimation  of  the  geostrophic  wind  that 
in  turn  leads  to  an  overestimation  of  the  upwell- 
ing  index.  The  assumption  made  in  presenting 
these  series  as  indicators  of  time  variations  of 
intensity  of  coastal  upwelling  is  that  the  actual 
pressure  gradient  at  the  coast  varies  in  rough 
proportion  to  the  total  onshore-offshore  pressure 
difference  which  is  used  in  the  computations. 

MEAN  YEARLY  CYCLE  OF 
INDICATED  UPWELLING 

Having  been  designed  to  give  a  consistant  in- 
dication of  temporal  variations  at  each  particular 


LOW 
PRESSURE 


Figure  5.  — Height  of  a  constant  pressure  surface.  The  slope  of  the  dashed  line  indi- 
cates the  assumed  gradient  used  in  the  calculations.  The  slope  of  the  solid  line  indicates 
the  lesser  gradient  which  mav  exist  in  realitv. 
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Table  3.  — Mean  monihly  values  of  the  indices  for  the  20-yr  period,  1948-1967.  Units  are  cu- 
bic meters  per  second  per  100  m  of  coastline. 


Jan. 

Feb. 

Mar. 

Apr. 

May 

June 

July 

Aug. 

Sept. 

Oct. 

Nov. 

Dec. 

60°N.149°W 

-138 

-  86 

-46 

-11 

0 

6 

6 

6 

-  3 

-26 

-73 

-109 

60°N,146°W 

-180 

-103 

-48 

-12 

.  2 

6 

5 

3 

-  9 

-34 

-94 

-129 

57°N,137°W 

-212 

-117 

-51 

-24 

-11 

0 

1 

-6 

-29 

-88 

-140 

-163 

54°N.134°W 

-  97 

■  68 

-27 

-20 

-10 

1 

3 

-1 

-V, 

-82 

-  98 

-  91 

51°N,131°W 

-  64 

-  36 

-12 

-  5 

4 

13 

16 

12 

-  3 

-40 

-  58 

-  57 

48°N,125°W 

-  90 

-  47 

-21 

0 

18 

25 

34 

22 

A 

-39 

-  88 

-100 

45°N.125°W 

-  94 

-  47 

-15 

9 

34 

48 

74 

50 

16 

-20 

-  73 

-  93 

42°N,125°W 

-  67 

-  28 

3 

33 

79 

103 

132 

91 

36 

0 

-  42 

-  57 

39°N,!25°W 

-  13 

9 

36 

69 

124 

168 

182 

139 

63 

20 

-     7 

-  12 

36°N,122°W 

11 

35 

80 

121 

203 

239 

198 

183 

94 

49 

12 

1 

33°N,1I9°W 

19 

48 

120 

178 

282 

312 

231 

212 

137 

76 

22 

10 

30°N,119°W 

56 

77 

116 

1-41 

199 

199 

143 

142 

129 

103 

65 

54 

27°N.116°W 

71 

93 

119 

148 

202 

195 

114 

105 

110 

106 

74 

63 

24°N.113°W 

51 

74 

93 

116 

143 

129 

48 

44 

49 

69 

52 

39 

21°N.!07°W 

18 

39 

97 

KMJ 

87 

39 

3 

5 

-14 

-15 

8 

8 

locution,  this  set  of  indices  may  be  admirably 
suited  for  describing  the  spatial  variations  of 
timing  and  duration  of  upwelling.  In  Table  3  are 
listed  the  long-term  monthly  mean  values  of  the 
indices  for  the  20-yr  period,  1948-67.  Figure  6 
displays  plots  of  these  values. 

The  period  of  maximum  values  of  the  indices  is 
earliest  to  the  south,  becoming  progressively 
later  in  the  year  with  distance  to  the  north.  At 
lat.  21  °N  the  peak  is  in  April.  Off  the  coast  of 
Baja  California  the  peak  is  in  May.  A  June  peak 
occurs  off  southern  and  central  California.  North 
of  San  Francisco  the  peak  is  generally  in  July, 
with  a  slight  shift  back  to  June  at  the  northern 
extremity  of  the  Gulf  of  Alaska. 

From  the  southern  tip  of  Baja  California  north- 
ward to  Monterey  Bay  the  mean  values  of  the 
indices  indicate  upwelling  throughout  the  year. 
Northward  from  San  Francisco  the  season  of 
indicated  upwelling  becomes  progressively  re- 
stricted, lasting  March  through  October  off  north- 
ern California,  April  through  September  off  Ore- 
gon, May  through  August  off  northern  Vancouver 
Island,  finally  narrowing  at  lat.  57 °N  to  a  scarce- 
ly significant  positive  value  of  the  index  during 
July. 

Tabulations  of  the  deviations  from  these  long- 
term  mean  values  accompany  the  tabulations  of 
the  values  of  the  indices  in  the  Appendix.  Graph- 
ical plots  of  both  the  anomalies  and  the  values 
of  the  indices  are  included. 


DISCUSSION 

Detailed  analysis  and  interpretation  of  the 
time-.series  plots  contained  in  the  Appendix  is 
beyond  the  scope  of  this  report.  However,  it  may 
be  useful  to  point  out  several  examples  and  to 
speculate  briefly  upon  possible  consequences. 

Accompanying  the  plots  of  the  monthly  indices 
for  each  location  are  tables  of  quarterly  and  an- 
nual averages.  The  following  three  examples  are 
each  illustrative  of  phenomena  occurring  basi- 
cally on  one  of  these  three  different  time-scales. 
As  it  happens,  the  first  two  concern  negative 
values  of  the  index.  Negative  values  of  the  index 
indicate  accumulation  of  wind-transported  sur- 
face waters  at  the  coast  and  resulting  down- 
welling. 

Short-Term  Intense  Convergences  in 
the  Northern  Gulf  of  Alaska 

A  glance  at  the  first  chart  in  the  Appendix 
reveals  an  extreme  negative  spike  in  the  up- 
welling index  at  lat.  60 °N,  long.  149  °W  during 
February  1950.  In  this  case,  the  indicated  con- 
vergence is  over  4y2  times  the  20-yr  mean  value 
for  February.  It  is  interesting  to  speculate  upon 
the  effect  of  such  an  anomaly  on  the  density  struc- 
ture of  the  northern  portion  of  the  Alaskan  Gyre. 

If  the  given  numerical  value  is  accepted  as  an 
estimate,  we  find  for  February  1950  an  anomalous 
onshore  flow  of  nearly  300  m-'  /sec  along  each 


Fisure  6. -Mean  mmithly  values  of  the  computed  upwelling  indices  for  the  20-yr  period,  1948  through 

through  1967  Im'"    sec'Vper  100  m). 
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100  m  of  coastline.  Assuming  this  flow  accumu- 
lates within  a  zone  extending  300  km  offshore, 
this  represents  for  the  1-mo  period  a  lens  of  ex- 
cess surface  water  about  24  m  thick.  If  the  flow 
away  from  the  coast  which  must  balance  this  ac- 
cumulation occurs  below  the  pycnocline  the  re- 
sult is  a  downward  shift  of  the  upper  density 
structure  an  average  of  24  m  over  the  shoreward 
300  km.  Such  a  shift  implies  an  extreme  increase 
in  baroclinicity  in  this  already  highly  baroclinic 
area,  tending  to  accelerate  the  westward  geo- 
strophic  flow  along  the  northern  boundary  of  the 
Gulf  of  Alaska. 

The  striking  negative  anomalies  for  the  same 
month  in  charts  of  the  index  at  lat.  60 °N,  long. 
146 °W  and  at  lat.  57  °N,  long.  137 °W  indicate 
that  this  feature  existed  over  the  major  portion 
of  the  northern  Gulf  of  Alaska  coastline.  This 
wide  geographical  extent  justifies  such  a  two- 
dimensional  treatment  as  presented  in  the  pre- 
vious paragraph.  It  also  indicates  a  possibly  far- 
reaching  effect  on  the  marine  environment  of  this 
region. 

Several  similar  short-term  intense  negative 
values  of  the  index  appear  in  the  charts  at  other 
times.  In  each  case  the  major  pulse  is  confined 
to  only  one  monthly  value  but  has  spatial  coher- 
ence in  the  Northern  Gulf.  Most  impressive  along 
with  February  1950  just  discussed  are  the  months 
of  January  during  the  years  1959,  1966,  and  1971. 
January  1971  differs  from  the  others  in  that  the 
absence  of  a  strong  anomaly  at  lat.  57°N,  long. 
135°W  indicates  a  lesser  eastward  extent  of 
the  feature. 

Periods  of  Intense  Winter  Convergence 
off  Washington  and  Oregon 

The  monthly  charts  of  the  index  at  lat.  45°N, 
long.  125 °W  indicate  the  periods  of  December 
1957  through  February  1958  and  January  through 
March  1961  as  having  been  intervals  of  extreme 
convergence  of  wind-drifted  surface  waters  at 
the  coast.  Winter  of  1958,  in  particular,  exhibits 
a  striking  anomaly  in  the  charts  for  locations 
extending  from  lat.  51  °N  to  lat.  39 °N.  These 
periods  differ  from  the  intense  convergences 
previously  described  for  the  more  northerly  re- 
gion in  containing  several  successive  months  of 
anomalously  large  negative  values  of  the  index 
rather  than  a  single  extreme  month   and   are 


highly  evident  in  the  charts  of  quarterly  sum- 
maries presented  in  the  appendix.  The  effect  of 
this  apparent  downwelling  at  the  coast  upon  the 
density  structure  would  be  to  accelerate  north- 
ward flow  or  to  decelerate  southward  flow.  Either 
situation  would  favor  an  anomalous  warm  ad- 
vection. 

Long-Term  Groupings  of 
Mean  Yearly  Values 

In  examining  these  data  for  long-term  or  re- 
gional trends  it  is  sometimes  helpful  to  filter  out 
some  of  the  shorter  term  variability  by  means 
of  annual  averages.  For  example  the  charts  of 
mean  yearly  values  indicate  anomalously  high 
mean  annual  upwelling  throughout  the  most 
recent  8-yr  period,  1964  through  1971,  from  Cape 
Blanco  to  San  Diego  (lat.  42  °N  to  lat.  33 °N).  In 
contrast,  along  Baja  California  (lat.  30  °N  to  lat. 
24  °N)  strikingly  low  anomalies  are  indicated  for 
the  4  successive  years,  1965  through  1968. 

The  period  1955  through  1959  has  anomalously 
high  values  from  the  Los  Angeles  Bight  to  the 
southernmost  point  covered  by  the  grid  (lat. 
33°N  to  lat.  21°N).  In  fact,  the  anomalously  high 
values  during  this  period  extend  all  the  way  north 
to  the  Straits  of  Juan  de  Fuca  (lat.  48°N)  if  the 
year  1958  which  is  anomalously  low  in  the  more 
northern  area  is  left  out.  The  low  yearly  mean 
values  for  1958  are  largely  due  to  the  extremely 
low  first  quarter  value  previously  discussed. 

Low  annual  absolute  values  of  the  indices  are 
found  for  the  first  4  yr  of  the  series  (1946-49)  at 
nearly  all  locations.  As  previously  mentioned  this 
may  be  due,  at  least  in  part,  to  relatively  sparse 
data  coverage  in  the  early  years  leading  to  some 
smoothing  of  the  gradients  used  in  the  computa- 
tions. 

Conclusion 

Anomalies  of  the  type  discussed  may  have  im- 
portant marine  biological  consequences.  Vertical 
displacements  of  the  thermocline  due  to  intense 
convergences  or  divergences  such  as  that  indi- 
cated for  February  1950  may  affect  the  depth 
on  the  continental  shelf  where  groundfish  are 
located.  Anomalous  upwellings  and  downwellings 
undoubtably  affect  productivity  and  the  concen- 
trations of  primary  producers  and  zooplankton. 
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They  may  also  affect  recruitment  and  distribu- 
tion of  fishery  resources,  particularly  where  the 
anomaly  is  of  longer  term  or  recurs  in  successive 
years. 

In  view  of  these  possibilities  it  is  hoped  that 
biologists  will  find  the  indices  presented  here  to 
be  convenient  and  suitable  for  incorporation  into 
studies  of  the  dynamics  of  biological  distributions, 
and  that  empirical  associations  may  lead  to  eco- 
logical hypotheses  and  finally  to  experiments 
which  can  verify  both  the  hypotheses  and  the 
indices  presented  in  this  report.  Finally,  although 
this  work  was  directed  primarily  toward  develop- 
ing a  tool  for  biological  and  fisheries  research,  it 
is  suggested  that  these  indices  may  find  appli- 
cations in  such  fields  as  coastal  oceanography, 
fog  research,  and  studies  of  climate  and  weather 
of  coastal  regions. 
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APPENDIX 


Tabulations  and  plots  of  coastal  upwelling  indices 


The  following  pages  display  the  coastal  upwell- 
ing indices  by  month,  quarter,  and  year  for  15 
near-coastal  points  over  the  period  1946-71.  The 
display  for  each  location  comprises  six  pages. 
The  first  four  contain  monthly  values.  The  fifth 
contains  quarterly  values.  The  sixth  contains  the 
completion  of  the  table  of  quarterly  values  and 
the  table  of  yearly  mean  values. 

The  first  column  of  figures  on  each  page  indi- 
cates the  month,  quarter,  or  year  to  which  the 
information  appearing  on  the  same  line  pertains. 
In  the  tabulations  by  month,  the  first  two  digits 
refer  to  the  year  and  the  last  two  to  the  month. 
For  example,  4601  refers  to  January  1946.  In 
the  tabulations  by  quarter  the  first  two  digits 
again  refer  to  the  year  and  the  last  two  refer  to 
the  quarter.  For  example,  4601  in  this  case  refers 
to  first  quarter  (Jan.,  Feb.,  Mar.)  1946.  In  the 
tabulations  by  year  the  year  is  printed  in  full. 

The  second  column  contains  the  value  of  index 
for  each  period.  The  units  are  metric  tons  per 
second  per  100  m  of  coastline.  These  units  may 
be  thought  of  as  indicative  of  the  average  amount 
(metric  tons  or  cubic  meters)  of  water  upwelled 
through  the  bottom  of  the  Ekman  layer  each  sec- 
ond along  each  100  m  of  a  straight  line  directed 


along  the  dominant  trend  of  the  coast  on  a  scale 
of  about  200  miles.  However,  it  must  be  remem- 
bered that  because  of  uncertainty  in  exact  values 
of  some  of  the  constants  employed  and  for  other 
reasons  outlined  in  the  section  on  "Method," 
these  indices  should  be  regarded  more  as  indic- 
ative of  relative  flucuations  of  upwelling  intensity 
at  each  point  than  as  quantitative  measures  of 
absolute  magnitude.  The  quarterly  and  yearly 
values  are  averages  of  the  monthly  values 
weighted  by  the  number  of  days  in  each  month. 
An  exception  to  this  rule  is  February  which  is 
considered  to  always  contain  28  days  in  order 
to  avoid  a  different  set  of  long-term  means  for 
leap  years.  Such  refinement  is  well  below  the 
level  of  sensitivity  of  these  indices. 

The  third  column  contains  the  anomalies  or 
deviations  from  the  long-term  mean  value  of  the 
index  at  the  particular  point  for  the  particular 
month,  quarter,  or  year  over  the  20-yr  period, 
1948-67. 

The  values  and  anomalies  are  plotted  to  the 
right  of  each  page.  The  heavy  line  indicates  the 
values  of  the  upwelling  index  while  the  lighter 
line  indicates  the  anomalies.  The  deviations  of 
the  anomaly  line  from  the  zero  line  are  shaded. 
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